Inversion of K3C60 Reflectance Data 



OO 
On 
G\ 

IX, 
(N 



(N 
O 
<N 
O 
OO 

o> 

-I— > 



o 
o 



X 



F. Marsiglio 1 ' 3 , T. Startseva 2 , and J. P. Carbotte 2 ' 3 
1 Department of Physics, University of Alberta, Edmonton, AB T6G 2J1 
2 Dept. of Physics & Astronomy, McMaster University, Hamilton, Ontario L8S J^Ml 
3 Canadian Institute for Advanced Research, McMaster University, Hamilton, ON L8S 4M1 

We outline a procedure for obtaining the electron-phonon spectral density by inversion of optical 
conductivity data, a process very similar in spirit to the McMillan- Rowell inversion of tunelling data. 
We assume both electron-impurity (elastic) and electron-phonon (inelastic) scattering processes. 
This procedure has the advantage that it can be utilized in the normal state. Furthermore, a very 
good qualitative result can be obtained explicitly, without iteration. We illustrate this technique on 
recently acquired far- infrared data in K^Ceo- We show that the electron-phonon interaction is most 
likely responsible for superconductivity in these materials. 



The determination of the underlying interactions and 
how they govern the various phases in materials (metallic, 
insulating, superconducting, etc.) is a major goal in cur- 
rent research of new materials. In the last ten years, both 
high-T c oxides and alkali-doped "bucky-balls" (A 3 C 6 o, 
etc.) have provided fascinating cases for study, since sev- 
eral phases are attainable by both doping and temper- 
ature variation. The general problem of understanding 
how these various phases come about is a formidable task, 
however, and remains unsolved. 

In this paper we investigate the origin of superconduc- 
tivity in K3C60, a superconducting fullerene with T c ss 19 
K. While several experimental facts point towards a con- 
ventional electron-phonon mechanism (eg. isotope effect 
measurements [0J|] and others — see Refs. ||, j|, j| and 
H for reviews), the key experiment, which is single elec- 
tron tunneling at biases above the gap edge, has not yet 
been performed to the extent that reproducible struc- 
ture is observed. In fact, definitive results from single 
electron tunneling may be very difficult to achieve for two 
reasons. First, the electron-phonon coupling strength 
may be very weak, in which case the structure will be 
difficult to observe, and second, neutron scattering exper- 
iments indicate that phonon modes are present at very 
high frequencies |8| (u> > 200 meV). If these modes are 
coupled to the itinerant electrons, then they may be im- 
possible to see through tunneling experiments, whose ac- 
curacy becomes less reliable at higher biases. 

For these reasons we attempt here to provide as defini- 
tive a demonstration of electron-phonon superconductiv- 
ity as tunneling normally provides, using, however, a 
higher frequency and potentially more sensitive probe, 
infrared spectroscopy. This idea is certainly not new, as 
many years ago Allen j^] provided a rudimentary basis for 
inferring the electron-phonon spectrum for lead, based on 
the absorption measurements of Joyce and Richards [10|. 
This work was followed up by Farnworth and Timusk [ 1 1 1 ; 
they utilized weak coupling expressions due to Allen and 
were able to extract a spectrum which agreed very well 
with that obtained from tunneling. 

The present work differs somewhat from these ap- 
proaches in that we utilize only normal state reflectance 



measurements, over the (in principle) entire frequency 
range. Through a Kramers-Kronig transformation we ob- 
tain the optical conductivity, from which we can deduce 
the electron-phonon spectrum. We first present some of 
the theoretical aspects, carry out some checks on Pb, 
and finally apply our procedure to K^Ceo- Our conclu- 
sion is that the optical measurements of Degiorgi et al. 
]l2[ are consistent with electron-phonon coupling of suf- 
ficient strength to drive the superconducting transition 
in K 3 C 60 . 

The starting point is the Kubo formula for the optical 
conductivity pq ] of electrons in a normal metal, at zero 
temperature: 
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is the electron self-energy due to the electron-phonon in- 
teraction (the term proportional to the electron-impurity 
scattering rate is explicitly included in Eq. (|l|)). Eq. (|l|) 
has been derived using a number of well-documented ap- 
proximations: impurity scattering is included in the sim- 
plest Born approximation, and electron-phonon vertex 
corrections have been ignored. Actually we have written 
the self-energy as a functional of a 2 F(il), although Allen 
Pj has noted (following Scher Jl4[ ) that this form can be 
used provided a 2 F(£l) is replaced by a more complicated 
spectral function called cq r F(SY), which in some limits is 
very well approximated by a 2 r F(fl). This latter function 
contains the essential "1 — cos(6)" factor for transport 
properties. So, in fact the above equations should con- 
tain the transport spectral function, but since this and 
the usual spectral function are very often qualitatively 
similar, we have dropped the distinction. 

It should be clear from Eqs. (|l|) and (|^) that the nor- 
mal state conductivity contains an image of the underly- 
ing electron-phonon spectral function. This is most ap- 
parent if we invoke a weak coupling approximation and 
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expand the denominator in Eq. ([j]) to first order in the 
self-energy, so one of the integrals can be performed by 
parts (l6) . Further straightforward manipulation leads to 
the following explicit expression for the electron-phonon 
spectral function: 
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This equation makes it clear that the extraction of 
a 2 F{v) by optical measurements is a difficult task. The 
data must be exceedingly accurate as two derivatives are 
required. Furthermore, the real and imaginary parts of 
the optical conductivity must themselves be obtained 
through Kramers-Kronig integrations of the reflectance 
(at least in the case of rcflcctometry) . This process re- 
quires extrapolations at low and high frequencies. Fi- 
nally, the plasma frequency must be obtained, either 
through a convenient sum rule fl5|| , or from an indepen- 
dent measurement. In practice we anticipate that accu- 
racy of reflectance data will limit this inversion technique 
to qualitative features of the spectral function. This lat- 
ter observation further justifies some of the approxima- 
tions made above concerning the neglect of vertex cor- 
rections, as they are anticipated to affect the result only 
quantitatively. 

It is instructive to determine the qualitative accuracy 
of Eq. (H) by comparing the result obtained with the 
true underlying spectral function. To this end we have 
carried out the following theoretical exercise: the con- 
ductivity for Pb has been evaluated using Eq. ([!]), start- 
ing with the McMillan-Rowell tunneling density of states, 
and then Eq. (||) is utilized to provide a 2 F(v). We have 
used several impurity scattering rates to illustrate that 
the result is largely independent of impurity scattering 
rate, as Eq. (|3|) implies. These results are plotted in 
Fig. 1. Also plotted are two curves which are essen- 
tially indistinguishable: one is the a 2 F(v) used in the 
calculation, the other is obtained through a numerical 
inversion of Eqs. (0) and (|J), without making the weak 
coupling approximations that led to Eq. (||) . The details 
of this more rigorous procedure are described elsewhere 
jl7j . We include the result here to show that the proce- 
dure works, although the problem is ill-defined (because 
of the need to unravel two integrations in Eq. ([j]) and 
(0)). We also remark that the accuracy requirements for 
the data are quite stringent, so that in the case of real 
reflectance data numerical inversion is problematic. Fig. 
1 illustrates that the explicit expression (^) is an excel- 
lent approximation to the true electron-phonon spectral 
function, particularly since Pb is a fairly strongly cou- 
pled metal (A = 1.5). We furthermore learn that using 
Eq. (^) will lead to unphysical negative components of 
a 2 F(v) at high frequency, but that these should essen- 
tially be ignored (the numerically inverted result does not 
contain these negative pieces, nor the spurious positive 
tails). It is also instructive to examine the effect 
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FIG. 1. The a 2 F{v) for Pb (solid curve) plotted vs. 
frequency, along with the estimates obtained from Eq. 
(0) with an impurity scattering rate, 1/r = (a) 1 meV 
(dotted), (b) 10 meV (dot-dashed), and (c) 20 meV 
(short-dashed-long-dashed). Note that all are qualitatively 
quite accurate, before they become negative at higher fre- 
quencies. Also plotted is the result (dashed curve, barely dis- 
cernible from the solid curve) obtained from a full numerical 
inversion, as described in the text. 
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FIG. 2. The o?F(v) for Pb (solid curve) plotted vs. fre- 
quency, along with the estimates obtained from Eq. (^|), us- 
ing 1/r = 10 meV, calculated at (a) T = 1, (b) T = 3.5, (c) 
T — 7, and (d) T = 14 K. As the temperature is raised, the 
fine details are lost, but the overall energy scale remains. 



of finite temperature, especially since we are exam- 
ining a normal state property, and the normal state is 
often inaccessible at low temperatures (because of super- 
conductivity, for example). In Fig. 2 we show the result 
of applying Eq. (||) to the conductivity for Pb calcu- 
lated at a scries of finite temperatures. The result illus- 
trates that at sufficiently high temperature the structure 
is lost, although the qualitatively correct result remains: 
that structure in the conductivity is due to inelastic scat- 
tering in the phonon region (in the case of Pb). This 
procedure works fairly well, even though Eq. (^ does 
not even follow in the weak coupling limit at finite tem- 
perature. With these limitations in mind, we proceed for 
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the remainder of this paper using the explicit equation 
©■ 

We begin by using the reflectance data at 25 K for 
K3C60 taken by Degiorgi et al. jli^ ]. We have performed 
our own Kramers-Kronig analysis on this data, to extract 
the complex conductivity. We should caution the reader 
that this procedure requires low frequency and high fre- 
quency extrapolations. Our results are essentially insen- 
sitive to the form of the high frequency extrapolation. 
The standard procedure for the low frequency regime is 
to use the Hagen- Rubens form for the reflectance [ fL8| |; 
however, we have used the full Drude form, since second 
derivatives of the conductivity are required, and hence 
the result may be sensitive to these approximations. We 
find (see below) that the low frequency result in partic- 
ular is indeed sensitive to our choice of low frequency 
extrapolation. 

To proceed with a low frequency Drude extrapolation, 
two parameters are required, the dc resistivity (as in the 
Hagen- Rubens form), and an effective electron scatter- 
ing rate, which we take to be independent of frequency 
for low frequencies. Once these are chosen (by matching 
the reflectance with the experimental result at some low 
frequency), the frequency-dependent phase can by deter- 
mined by a Kramers-Kronig integral . From these fol- 
low all the optical properties of the solid, in particular the 
real and imaginary parts of the conductivity. Equation 
(|J) then provides us with an estimate for a 2 F(v). The 
result is plotted in Fig. 3, along with the neutron scatter- 
ing data from We note several features in the figure. 
First, there are indeed high frequency negative regions, 
which our experience with Pb teaches us to ignore. Sec- 
ond, there is a prominent low frequency negative region 
(near 20 meV) for which we have no explanation, except 
that perhaps the low frequency extrapolation may be in 
error. Also note the very low frequency peak at about 
5 meV, which is actually in agreement with the neutron 
data. On the other hand, an analysis using the Hagen- 
Rubens low frequency form (which is less accurate than 
the Drude form) leads to essentially the same result mi- 
nus the low frequency peak. Hence, we cannot say for 
sure that this peak (attributed to librons) is coupled to 
the electrons. The negative region just above the libron 
peak persists regardless of the method of low frequency 
extrapolation; aside from the point noted above, it could 
be due to varying electron density of states effects [[l9| , 
for example. We proceed, assuming the problem can be 
resolved by including such corrections, and discard the 
negative pieces in Fig. 3 for the remainder of the paper. 

If we ignore the negative pieces of a 2 F(v) the quali- 
tative agreement with the neutron scattering data (plot- 
ted with a dashed curve [|| ) is striking. Not only is 
the energy scale correct, but peaks line up fairly well, 
suggesting that the coupling is more or less energy in- 
dependent. It is difficult to assess whether contributions 
are truly coming from the higher energy region, as the 



neutron scattering data ends around 200 mev. 

In any event, we can think of the a 2 F(y) in Fig. 3 
as an " educated guess" and now proceed to calculate the 
expected reflectance. One might first try to fit the data 
by selecting a suitable plasma frequency and frequency- 
independent (impurity) scattering rate, so that the con- 
ductivity is Drude-like. Such an attempt fails, and it is 
clear that inelastic scattering is required. We use Fig. 3 
without adjustment and vary (by hand) the plasma fre- 
quency and impurity scattering rate to obtain reasonable 
agreement with the measured reflectance. The result is 
plotted in Fig. 4. While significant changes in a 2 F{y) 
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FIG. 3. The a 2 F(v) for K 3 C 6 o (solid curve) extracted 
from the reflectance data of Degiorgi et al. using Eq. 

(^). For purposes of analysis we have omitted the negative 
parts (see text). Also shown are the neutron scattering re- 
sults from Ref. || (dashed curve) . Clearly the energy scale in 
a 2 F{v) matches that of the phonons, and some of the peaks 
even line up correctly. 
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FIG. 4. The reflectance calculated using the spectral func- 
tion from Fig. 3 (solid curve) along with the drude parameters 
up — 6000 cm -1 and 1/r = 95 meV. The reflectance data of 
Degiorgi et al. Jl2| (dashed curve) are also shown. 

will lead to a deterioration of the agreement shown in 
Fig. 4, minor adjustments can be imposed which will ac- 
tually improve the agreement. However, this latter fine- 
tuning violates the spirit of Eq. (p|), which, as we learned 
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from the exercise with Pb, provides a qualitatively cor- 
rect estimate for the bosonic spectral function. 

We should note that we obtain as fitted parameters 
the result u>p = 6000 cm -1 and 1/r = 95 meV for the 
plasma frequency and impurity scattering rate, respec- 
tively. The first result is somewhat less than the result 
obtained via the sum rule (TJ] , which was used in Eq. (||) 
(ujp = 7200 cm -1 ). We view this as a rather minor in- 
consistency, compared to, say, that which exists between 
the penetration depth-derived plasma frequency, and the 
one from the optical sum rule jl2). The value obtained 
for the impurity scattering rate is a little more discon- 
certing, since it is quite large. In the approach of Rcf. 
|l2f such a large scattering rate is avoided by explicitly 
accounting for a mid-infrared band. We could, of course, 
also include a mid-infrared band with extra parameters 
to improve the overall fit and achieve a more realistic 
scattering rate, but then the derived a 2 F(v) would have 
more model dependence then is already implied in Eqs. 
©and®. 

Given the extracted a 2 F(y) in Fig. 3 (without the 
negative components) one can ask whether or not such 
a spectral function can account for the superconduct- 
ing properties in 'K^Cqq. To this end we have com- 
puted the Coulomb repulsion pseudopotential, /i*, given 
that T c = 19 K. The result is fi* = 0.4 (using a cut- 
off, uj c = I eV). The spectrum can be characterized by 
A = 1.2 and uj hl = 40 meV. This gives us T c /uji n = 0.04, 
which, in spite of the rather large value of A, represents 
a fairly weak coupling strength. Nonetheless, the gap ra- 
tio, 2Ao/fcsT c comes out fairly close to 4, significantly 
different from the BCS weak coupling result, a result 
which is partially due to the large value of n* p3). If 
the low frequency peak in a 2 F(v) at 5 meV is excluded, 
the results are then very BCS-like: A = 0.8, ^* = 0.34, 
T c /uji n w 0.016, and the gap ratio is 2A /fcsT c w 3.6. 
Either of these possibilities is strongly refuted by a gap 
measurement of order 5 [^lj , though it is difficult to dis- 
criminate between the two lower values based on far- 
infrared measurements. Microwave measurements may 
provide a discriminating probe, as discussed in p2fl. 

In conclusion, we have described an inversion scheme 
to extract the source of inelastic scattering in a metal. 
Such a scheme has many advantages over that already 
utilized in tunneling, not the least of which is Eq. (||), 
which allows a direct estimate from optical measure- 
ments with essentially no effort. Fairly accurate experi- 
ments are required however. We applied this technique to 
K3C60, and have found the qualitatively correct a 2 F(y) 
for this material. The coupling is sufficiently strong to 
explain the superconductivity in this material; in fact 
the n* required is conspicuously large, so that, as sus- 
pected by many, electron-electron correlations are strong 
as well. Nonetheless, the weak coupling approach pre- 



sented here appears to be able to account for the main 
features observed in the far-infrared. K3C60 is a con- 
ventional weak coupling electron-phonon superconductor 
with an anomalously large Coulomb pseudopotential. 

We acknowledge helpful discussions with Tom Timusk. 
This research was supported by the Natural Sciences 
and Engineering Research Council (NSERC) of Canada 
and by the Canadian Institute for Advanced Research 
(CIAR). 



[1] C.C. Chen and CM. Lieber, Science 259, 655 (1993). 
[2] A. P. Ramirez et al. Phys. Rev. Lett. 68, 1058 (1992). 
[3] M.P. Gelfand, Superconductivity Review 1, 103 (1994). 
[4] W.E. Pickett, Solid State Physics 48, 226 (1994). 
[5] C.H. Pennington and V.A. Stenger, Rev. Mod. Phys. 68, 
855 (1996). 

[6] O. Gunnarsson, Rev. Mod. Phys. 69, 575 (1997). 

[7] J. Ostrick has performed some preliminary tunneling 
measurements on these systems (private communica- 
tion) . 

[8] L. Pintschovius, Rep. Prog. Phys. 57, 473 (1996). 
[9] P.B. Allen, Phys. Rev. B 3, 305 (1971). 
[10] R.R. Joyce and P.L. Richards, Phys. Rev. Letts. 24, 1007 
(1970). 

[11] B. Farnworth and T. Timusk, Phys. Rev. BIO, 2799 
(1974); ibid B14, 5119 (1976). 

[12] L. Degiorgi, E.J. Nicol, O. Klein, G. Griiner, P. Wachter, 
S.-M. Huang, J. Wiley and R.B. Kaner, Phys. Rev. B49, 
7012 (1994). More recent results are presented in Degiorgi 
et al. Nature 369, 541 (1994). 

[13] G.D. Mahan, Many-Particle Physics (Plenum Press, New 
York, 1981). 

[14] H. Scher, Phys. Rev. Letts. 25, 759 (1970). 

[15] R. Kubo, J. Phys. Soc. Jpn. 12, 570 (1957). 

[16] We perform the weak coupling expansion somewhat dif- 
ferently from Allen || . The denominator v + £ is kept as 
the unperturbed part; once the integral in the expanded 
portion is performed, we de-expand (as Allen does — see 
his Eq. (61)) but we do not attempt to extract a renor- 
malization component which leads to his undetermined 
function $(u). Our result has no unknown components. 

[17] F. Marsiglio, unpublished. 

[18] F. Wooten, Optical Properties of Solids (Academic Press, 

New York, 1972). 
[19] B. Mitrovic and M.A. Fiorucci, Phys. Rev. B31, 2694 

(1985). 

[20] C.L. Sammer and J. P. Carbotte, Phys. Rev. B48, 3375 
(1993). 

[21] Z. Zhang, C. Chen and CM. Lieber, Science 254, 1619 
(1991). 

[22] F. Marsiglio, J. P. Carbotte, R. Akis, D. Achkir and M. 
Poirier, Phys. Rev. B50, 7203 (1994). See also F. Mar- 
siglio and J. P. Carbotte, Aust. J. Phys. 50, xxxx (1997). 



4 



